[x1,x2] = meshgrid(-2:0.1:3,-2:0.1:2);
 z1=10*(x2-x1.^2);
z2=1-x1;
z=z1.^2+z2.^2;
 mesh(x1,x2,z)
 xlabel('x1'),ylabel('x2'),zlabel('z')